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DNA-spools, structures in which DNA is wrapped and helically coiled onto itself or onto a protein 
core are ubiquitous in nature. We develop a general theory describing the non-equilibrium behavior 
of DNA-spools under linear tension. Two puzzling and seemingly unrelated recent experimental 
findings, the sudden quantized unwrapping of nucleosomes and that of DNA toroidal condensates 
under tension are theoretically explained and shown to be of the same origin. The study provides 
new insights into nucleosome and chromatin fiber stability and dynamics. 

PACS numbers: 87.15.He, 82.37.Rs, 87.16. Sr 



Introduction: Wrapped DNA-protein complexes are 
ubiquitous in nature (l| and play key roles in many fun- 
damental life processes. Prominent examples of DNA 
wrapping proteins are: the Lacl repressor [2$ participat- 
ing in the bacterial gene regulation, the DNA-gyrase J 
directing changes in DNA topology, RNA polymerase ^ 
copying DNA to RNA, and the histone octamer || per- 
forming DNA packaging into nucleosomes leading in each 
cell to the enormous condensation of meters of DNA into 
micron sized chromosomes. Besides the natural wrapped 
architectures there are attempts to design nanoparticles 
imitating that motive p as a means to efficiently pack 
and transport DNA into cells. In most of these ligand- 
DNA complexes the geometry and chemistry of the ligand 
surface enforces the DNA to follow a superhelical wrap- 
ping path with one or more tight turns. Remarkably, 
upon addition of multivalent condensing agents (like in 
sperm cells) or under high crowding conditions (like in 
virus capsids or during ^-condensation) DNA also shows 
an intrinsic ability to self-organize into large toroidal 
spools |7j. 

In the past decade single molecule experiments have 
become available allowing to apply tension to individual 
polymers in order to probe their mechanical prop erties || 
as well as their interaction with ligands [tj Hoi Hl| and 
molecular motors 01 • Static and dynamic force spec- 
troscopy |l3l| developed into a powerful tool for measur- 
ing equilibrium as well as kinetic characteristics of single 
molecules, going far beyond the information accessible by 
classical bulk experiments. Application of these methods 
to DNA-spool geometries has been awaited for l ong and 
was reported only recently for single nucleosomes [14( and 
single DNA toroidal condensates 0, 0] . These experi- 
ments - at first glance completely unrelated - reside on 
different length and energy scales and ground on differ- 
ent mechanisms of wrapping. Despite that, they both 
reveal the same surprising result apparently contradict- 
ing all the available bulk data: the unfolding of wrapped 
DNA from the spools is a catastrophic event, i.e., it is 
sudden and quantized and happens one DNA turn at a 
time. The aim of this letter is to theoretically explain 
this unusual nonequilibrium effect and to demonstrate 




FIG. 1: A) An unfolding DNA-protein spool under tension is 
characterized by two angles: the desorption angle a and the 
tilt angle /3. B) The energy landscape of a DNA-spool (as 
given by Eq.Q]for A = 50nm, F = ApN, e ads = 0.7k B T/nm, 
R = 4.2nm, H — 2.5nm) and its unfolding pathway. 

the universality behind it. Our theory is then applied 
to nucleosomes and DNA toroids allowing us to extract 
from experiments the relevant energetic parameters and 
to resolve apparent " oddities" in the dynamics of these 
systems. 

General Model: The DNA is assumed to be adsorbed 
on the protein spool surface along a predefined helical 
path with radius R and pitch height H. This path ac- 
counts for the typical chemical structure of such a pro- 
tein spool surface (e.g. distribution of charges, hydrogen 
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donors/acceptors etc.) along which the DNA is adsorbed 
with a net adsorption energy density e a d s given by the 
difference of the total DNA-protein binding energy and 
the stored DNA bending energy per length along the he- 
lical path. The degree of DNA adsorption is described 
by the desorption angle a which is defined to be zero 
for one full turn wrapped (cf. Fig. 1). After short in- 
spection it becomes clear that the unwrapping problem 
is non-planar and that the spool needs to rotate tran- 
siently out of the plane while performing a full turn - an 
effect already pointed out by Cui and Bustamante (lif . 
Therefore a second angle, /3, is introduced to describe the 
out-of-plane tilting of the spool. When a tension F (along 
the Y-axis) acts on the two outgoing DNA "arms" the 
system (i.e., the wrapped spool together with the free 
DNA ends) will simultaneously respond with (?) DNA 
deformation, with («) spool tilting and with (in) DNA 
desorption from the spool. In the following we assume 
that the DNA has freely rotating ends (as in the exper- 
iment [lj|) that allows us to neglect the twist degree of 
freedom. Then the total energy of the system as a func- 
tion of a and (3 writes E tot = 2Re a d s a + 2Eb en d — 2FAy. 
Here the first term is the adsorption energy, the second 
the bending deformation energy of the two free DNA por- 
tions, and the third term describes the gain in potential 
energy by pulling out the DNA ends by a distance 2Ay. 
Applying linear elasticity theory and elementary geom- 
etry E to t can be redistributed into the following three 
terms: 

Etot (a, 13) = E comp (a) + E geom (a, (3) + E sti ff (a, (3) 

(1) 

The first term E comp — 2R (e a ds — F)a describes 
the competition of the adsorption and the ap- 
plied force. The "geometrical" energy term 
E geom = 2Fi?(cos/3sina- [H / 2n R) (tt - a) sin f3) 
stems from the gain/loss of potential energy by spool 
opening (change of a) and rotation (change of (3). Fi- 
nally, the last and most remarkable term E st i ff= 8\/ AF 

(l - + (R/R) cos P cos a + (H/2ttR) sin/3) /2J 

accounts for the stiffness of the non-adsorbed DNA 
portions. Here A w 50fcsTnm is the DNA bending 
stiffness at room temperature and R = R 2 + H 2 /An 2 . 
Two effects contribute equally to E^iff' 1) the bending 
energy of the deformed DNA arms and 2) the loss 
of potential energy by "wasting" length due to DNA 
deformation. To understand the implications of Eq. ^ 
on the kinetics of unwrapping we consider two limiting 
cases. First let us look at the case of a large thin spool, 
i.e., R A/ksT (or, equivalently, an infinitely flexible 
polymer) and R 3> H, where we may neglect E st iff. 
In that case and for F > e a ds the spool moves from 
the (thermodynamically) metastable state Mi with 
a = a n = — arccos(l — e a ds/F) and (3 — via a saddle 
point S at a = and j3 = «o into a more favorable 



minimum M2 at a = ir + ctQ and (3 = 7r. Remarkably 
S constitutes a significant energetic barrier between 
Mi and M 2 given by AE to t = 2FR(a cosao — sinao). 
For hypothetical yet reasonable parameter values, say 
R = 50nm, e a d s — lksT /nm and F = 2e a ds we obtain 
a huge barrier of AE to t ~ lOksTl A second interesting 
limit of Eq. ^ is given by a flat spool and high polymer 
stiffness, i.e., A ^> Rkj^T and R^> H . For not too large 
forces (F < A/R 2 ) and e a ds S E the kinetic behavior 
is roughly dominated by the term E st iff. In this case 
we find a transition path from (a, [3) = (0, 0) over the 
saddle point (ir/2,ir/2) to the state (tt,it) with a barrier 
height AEtot = 8%/AF(l - 1/V2). Note that in this 
limit the DNA actively participates in the suppression 
of unwrapping (AE to t ~ A 1 / 2 F 1 / 2 ) which can even give 
rise to negative resistance effects 0] for small forces. 
In preliminary conclusion, in both limiting cases the 
unwrapping meets significant kinetic barriers but for 
different reasons: because of unfavorable projection of 
the force in terms of the (a, (3) configurational space in 
the first limit and due to significant transient bending of 
the DNA arms during the transition in the second limit. 
For realistic DNA-spools we are somewhere in between 
these two cases. 

Nucleosome Unwrapping: The most abundant DNA 
spool in nature is the nucleosome where 1 and 3/4 turns 
of DNA, 147 bp, are wrapped around a protein core on 
a lefthanded superhelical path with diameter 4.2nm and 
2.5nm pitch. The question about the equilibrium and 
kinetic stability of nucleosomes is one of the important 
experimentally unsettled questions in present molecular 
biology. How can nucleosomes be highly stable with 
its wrapped DNA being highly accessible at the same 
time 17]? A recently performed experiment |l4| mea- 
suring the critical force required to unwrap single nu- 
cleosomes reveals an interesting and unexpected behav- 
ior When small forces (F < lOpN) are applied for 
short times (~ 1 — 10 s) the nucleosome unwraps only par- 
tially by releasing the outer 60-70 bp of wrapped DNA 
(moving from state a to b in Fig. IB) in a gradual and 
equilibrium fashion. For higher forces (F ^ 20pN) nu- 
cleosomes show a pronounced sudden non-equilibrium re- 
lease behavior of the remaining 80 bp (cf. c-g in Fig. IB) 
- the latter force being much larger than expected from 
equilibrium arguments [l9j . In fact, experiments (l7]| 
measuring spontaneous partial unwrapping of nucleoso- 
mal DNA suggest 30/csT per 147bp leading to an unpeel- 
ing force of ~ 2.5pN. To explain this peculiar finding 
Brower-Toland et al. 0] conjectured that there must be 
a barrier in the adsorption energy located after the first 
70-80 bp which reflects some biochemical specificity of 
the nucleosome structure at that position. Their analy- 
sis of the dynamical force spectroscopy measurements re- 
vealed an apparent barrier of ~ 38/csT smeared out over 
not more than 10 bp. However, there is no experimental 
indication of such a huge specific barrier - neither from 
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B) Force F [pN] 

FIG. 2: Kinetic barriers opposing DNA-spool unfolding as 
function of applied tension computed for: A) The nucleo- 
some (R = 4.2nm, H = 2.4nra, s° ds = 0.7k B T/nm, A = 
50fc_BTnm,cf. text) for various interstrand repulsion energy 
densities e es and B) the DNA toroid (R = 50nm, H — 2Anm, 
A — 40fcsTnm) for various adsorption energy densities e a ds- 

the crystal structure [j| nor from the equilibrium accessi- 
bility to nucleosomal DNA ^t}- Consequently the ques- 
tion arises if the barrier is really caused by biochemistry 
of the nucleosome or, as we show below, by its underlying 
geometry and physics. 

To see that the effect is mainly physical we apply 
Eq.^to compute the barrier. For this purpose we model 
the nucleosomal adsorption energy density as e a ds (a) — 
e ads + @( a )£ es where e ads « 0.7fcsT/nrn is taken from 
the reversible part (for the first 60-70bp) of the measure- 
ment in Ref. |l4| . The introduction of the step function 
(6 = for a < and 8 = 1 for a > 0) together with 
a new parameter e es , the electrostatic interaction energy 
density, accounts for the DNA-DNA repulsion of the two 
adjacent helical gyres which acts only for a < reduc- 
ing the net e a ds (cf. Fig. 1A). Using e ads from above 
we can compute the barrier height for nucleosome un- 
folding for various values of e es as done in Fig. 2A. To 
relate the barrier heights from Fig. 2A to the dynami- 
cal force spectroscopy (DFS) measurements in Ref. 0] 
we generalize the classical relation between the loading 
rate rp and the most probable rupture force F* to 
the case of nonlinear force-barrier dependence and obtain 




FIG. 3: Optimal fits of the DFS data from Ref. for var- 
ious attempt frequencies give the corresponding electrostatic 
DNA-DNA-repulsion e es . 



In (r F /r ) = In [-u aU k B T/ {r AE')}- AE/k B T. Here r F 
and AE are functions of F* and AE 1 = d (AE) /OF*. 
ro = lpNs^ 1 is an arbitrary scale on the r F axis and 
Vatt is the typical attempt frequency of the nucleosome. 
Assuming v at t in the range 10 5 to 10 8 s _1 we can fit the 
experimental data from |l4| to obtain the correspond- 
ing values of e es , cf. Fig. 3. Keeping in mind that v aU 
is dominated by the slowest process involved in the un- 
folding event we estimate v att ~ 10 6 s _1 j2jj. The latter 
implies e es ~ 1.4— X.Dz^Tjnm. So indeed, at the line 
a = there is clear jump in adsorption energy density 
as we would naively expect from repulsive DNA-DNA 
electrostatics under these conditions 21]. This explains 
why under equilibrium conditions (at F — 0) the DNA 
deeply inside the nucleosomes (almost the whole bound 
DNA!) can be rather easily accessed by proteins [l7j but 
the nucleosome is still highly stable: The line a = 
can be moved to each position inside the nucleosome if 
the left and right DNA arms are adsorbed/desorbed in 
a consistent manner. Beyond that (a > 0) the assisting 
electrostatics switches off and the nucleosome is suddenly 
strongly stabilized (by 60 — 70/cbT in total!). 

Toroid Unspooling: When long DNA molecules con- 
densed with multivalent counterions were stretched in a 
single molecule experiment 10] Baumann et al. found 
a surprising behavior. When a critical force (typically 
F w 4 — 12pN) is reached large portions of DNA are 
released in packets in a discontinuous manner ("stick 
release pattern" ) . When the same experiment was re- 
done recently by another group [ill ] a pronounced quan- 
tization in the DNA release length of ps 300nm was 
clearly demonstrated. It was noted in Refs. [Tot Hlf 
that the latter correlates exactly with the typical size 
(R w 50nm) of toroidal condensates formed in solution 
and led those authors to the conclusion that a single 
turn of DNA unwraps from the toroid spool at a time. 
Despite that interesting finding the mechanism behind 
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this non-equilibrium effect remained unexplained. How- 
ever in the light of our theory the explanation is again 
straightforward as a DNA-toroid exhibits a spool geom- 
etry with R « 50nm and H -C R. The "first limit" 
considered above gives here a good approximation. The 
barrier heights for different values of e a ds as a function 
of force are computed in Fig. 2B. Similar as in the case 
of the nucleosome the attempt frequency v aU is domi- 
nated by the rotational friction - here of the 50nm sized 
toroid object - leading to v att « 3 x (10 2 — 10 3 ) s _1 . For 
high concentrations of the condensing a gen t spermidine 
one finds e ads w 0.2 - 0M B T/nm (cf. [3 Ell and the 
references therein). In case of equilibrium this means 
a very small peeling-off force of F sa 1 — 1.5pN. Our 
model together with Fig. 2B allows us now to predict the 
activated non-equilibrium behavior to have very low un- 
folding frequencies v un f — v att exp (-AE (F) /fcsT), for 
instance 1CT 6 to lCT 2 ^ 1 for F = 2pN , 10" 3 to Is" 1 for 
F = ApN and 0.3 to 50 s" 1 for F = 8pN ^ - consistent 
with experimental findings [Tollllj|. 



Conclusion: We have shown that DNA-spools rang- 
ing from protein-DNA complexes to DNA toroids share 
a universal feature inherited by their geometry: They 
are strongly kinetically protected from mechanical dis- 
ruption upon applied tension. In the case of chromatin 
fibers consisting of large arrays of nucleosomes and other 
DNA spooling proteins this effect provides a great biolog- 
ical advantage. Strong molecular motors like RNA poly- 
merase and helicase or microtubuli during cell division 
are known to act on the fiber with significant transient 
tensions of the order of 20pN or even more. While a hy- 
pothetical "fiber A" consisting of DNA and non-spooling 
proteins (say only DNA bending proteins) would imme- 
diately lose most of its protein content a "fiber B" con- 
stituted of DNA-spools would survive long time periods 
(up to 10 6 - 10 8 times longer than "fiber A"). We can 
speculate that this obvious advantage was not overlooked 
by nature and has flown into the chromatin fiber design 
and the nucleosome-spool shape. The remarkable uni- 
versality of the "kinetic protection" also shows up in the 
case of DNA-toroids which are roughly ~ 10 times larger 
while the DNA is ~ 10 times weaker adsorbed than for 
typical DNA-protein spools. While the biological impli- 
cations of this finding still have to be fully figured out it 
seems that this might play a role in the injection / ejec- 
tion process of DNA from viral capsids inducing similar 
quantization effects as found here. Looking at the wealth 
of peculiar effects revealed by the single molecule exper- 
iments 0, El we feel that the present understanding of 
DNA condensation kinetics is still incomplete, yet one 
partial mystery seems resolved. 
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sphere is of the order 2 
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- 10 f, s" 1 with rj s 
The typical fre- 



being the water viscosity (a centipoise) 
quency that characterizes the relaxation of the DNA arms 
is comparable to that - even if one accounts for addi- 
tional complexed nuclesomes as it is the case in the ex- 
periment p4l |. There the first unfolding nucleosome is sur- 
rounded by 16 other nucleosomes that have to move via 
a distance As ~ 25nm (unfolding length) under a force 
of F « 10-20piV leading to an lower bound 10 5 - 10 6 s _1 
of that frequency. 

Note that the latter alone would not explain the sudden 
catastrophic behaviour and the slow kinetics of unfolding 
if the effects described by Eq^were not included explic- 
itly. 

The upper/lower estimate correspond to e a ds = 0.2 / 
0.3fc_e T/nm; note the strong dependence of u un f on e a ds- 
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